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ABSTRACT 

We present AGN from the Sloan Digital Sky Survey (SDSS) having double- 
peaked profiles of [O ill] AA 5007, 4959 and other narrow emission-lines, mo- 
tivated by the prospect of finding candidate binary AGN. These objects were 
identified by means of a visual examination of 21,592 quasars at z < 0.7 in SDSS 
Data Release 7 (DR7). Of the spectra with adequate signal-to- noise, 148 spec- 
tra exhibit a double-peaked [O ill] profile. Of these, 86 are Type 1 AGN and 
62 are Type 2 AGN. Only two give the appearance of possibly being optically 
resolved double AGN in the SDSS images, but many show close companions or 
signs of recent interaction. Radio-detected quasars are three times more likely to 
exhibit a double-peaked [O ill] profile than quasars with no detected radio flux, 
suggesting a role for jet interactions in producing the double-peaked profiles. Of 
the 66 broad line (Type 1) AGN that are undetected in the FIRST survey, 0.9% 
show double peaked [O ill] profiles. We discuss statistical tests of the nature 
of the double-peaked objects. Further study is needed to determine which of 
them are binary AGN rather than disturbed narrow line regions, and how many 
additional binaries may remain undetected because of insufficient line-of-sight 
velocity splitting. 

Previous studies indicate that 0.1% of SDSS quasars are spatially resolved 
binaries, with typical spacings of ~ 10 to 100 kpc. If a substantial fraction of the 
double-peaked objects are indeed binaries, then our results imply that binaries 
occur more frequently at smaller separations (< 10 kpc). This suggests that 
simultaneous fueling of both black holes is more common as the binary orbit 
decays through these spacings. 
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Subject headings: galaxies: active — quasars: general — black hole physics 



Introduction 



Binary quasars are a rare but importa nt aspect of galactic evolution and the AGN phe- 
nomenon, as reviewed by iKomossal ( 120061 ). Studies of the nature and incidence of binary 
quasars provide insight into galaxy mergers at substantial look-back times and into the fuel- 
ing of AGN by accretion onto supermassive black holes. According to current cosmological 
models, many large galaxies have undergone at least one major merger. A likely mechanism 
for AGN fueling is the migration of gas caused by tidal torques in a merger. Therefore, one 
might expect to see a substantial number of binary quasars. 

Evidence for binary AGN includes peculiar morphologies of radio galaxies, X-ray re- 
solved double nuclei, and opticall y resolved double QS Os. Out of 10 5 QSOs, only ~ 10 2 
optical binary quasars are known ( [Foreman et al.ll2009l and references therein). Most have 
spacings of 10 kpc or more, and lie at z ~ 2 to z ~ 3, where study of the host galaxy is 
difficult. Binaries with smaller s eparations include a 2 kp c (0.3 arcsec) optical example at 
z = 0.848 observed with HST bv I Junkkarinen et al.l (boqiL a nd the ULIRG NGC 6240 with 
a double X-ray nucleus separated by 1.4 kpc ( IKomossal 120031 ) . Both authors argue that the 
observed incidence of such close binaries is roughly 100 times less than might be expected if 
the fueling of one black hole giving QSO activity at the relevant stages of a merger mandates 
the fueling of the second black hole. A single, m uch closer (7 pc) binary AGN candidate is the 
radio galaxy 0402+379 ( iRodriguez et al.l 120091 ) . Additional examples of optically resolved 
binary AGN include a z =0.36 galaxy from the COSMOS survey with a spacin g of 2.5 kpc 



(IComerford et al.ll2009bl ) and a z=0.44 binary quasar with a spacing of 21 kpc (IGreen et al. 
2010h. 



The narrow emission lines of AGN offer a potential indicator of binary AGN at inter- 
mediate separations of ~ 1 to 10 kpc. The radius of the narrow emission-line region (NLR) 
is typically a few hundred parsecs, so that such objects should have distinct NLRs orbiting 
with each of the merging nuclei with its black hole. Typical orbital velocities in galaxies are 
of order a few hundred kilometers per second, similar to the widths of AGN narrow lines. If 
two adjacent narrow-line regions in a merging system were detected by a single spectrograph 
slit or fiber, a double-peaked emission-line profile could result. This provides a technique for 
identifying candidate binary AGN th at are spatially unresolved. Several e xamples have al- 



ready b een discussed in the literature ( iZhou et al.ll2004l : iGerke et al.ll2007l ). IComerford et al. 



( I2009al ) have studied AGN with narrow lines displaced in velocity from the host galaxy as 
possible examples of galactic mergers in progress, including two objects with double-peaked 
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narrow lines. 



We have carried out a search for double-peaked narrow line profiles, primarily the [O ill] 
lines, in quasar spectra from the Sloan Digital Sky Survejl^]. The goals of this work were 
to identify interesting individual objects for study, and to constrain the frequency of binary 
quasars in the relevant regime of separation. This search has identified 148 double-peaked 
[O ill] objects, with velocity splittings between 180 km s _1 and 1400 km s" 1 . This includes 
roughly equal numbers of AGN 1 (broad lines) and AGN 2 (narrow lines only). Here we 
describe the search and the double-peaked objects, consider alternative causes of the double- 
peaked profiles, including bipolar jets and disk rotation, and discuss statistical inferences for 
AGN fueling in merging galactic nuclei. We assume a concordance cosmology with Ho = 70 
km s^Mpc" 1 , n A = 0.7. 



Recently, two independent studies by iLiu et al.l ( 120091 ) and IWang et all ( 120091 ) have 
appeared, involving narrow line AGN 2 selected from the SDSS galaxy database. These 
authors find a similar incidence of double [O ill] objects to that found here. The two studies 
are complementary to ours, in that our work emphasizes broad line objects. 



2. Sample and Method 

Our spectra are obtained from SDSS DR7, which contains 21,592 QSO spectra within 
our redshift range, 0.1 < z < 0.7. Because it is possible to have strong [O ill] emission 
despite a noisy spectrum, we imposed no a priori signal-to- noise (S/N) cutoff. We conducted 
a visual inspection of the spectra in the region of several key emission lines. Our primary 
search criterion was a double-peaked [O ill] profile in both A5007 and A4959, consistent with 
the 3:1 intensity ratio fixed by atomic physics. Once these objects had been identified, a 
more thorough inspection of each double [O ill] object was conducted with the IRAF routine 
SPLOT | to determine whether the velocity splitting observed in [O ill] was detectable in 
other emission lines. The lines considered were H/3 , [O ill] , [O n] , [Ne ill] , and [S n]. 
While confirmation of the [O ill] split in these lines strengthens the case, we did not require 
it for inclusion in the list of double-peaked objects if the S/N was such that the double-peak 
was not ruled out. For approximately 12 of our objects, the [Ne v] line showed sufficient 
S/N to allow measurement of the split. In these cases, there was agreement with the [O ill] 



1 The SDSS website is http://www.sdss.org 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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velocity split within about 10%, similar to the measurement accuracy. In a few more cases, 
the S/N allowed qualitative confirmation of the split by visual inspection. It appears that 
the split is seen in [Ne v] whenever the S/N is adequate to see it, Thus, the double-peaked 
profile extends to higher ionization lines in most if not all cases. However, [Ne v] was rarely 
useful to confirm an [O ill] split that was not already confirmed by other lines. 

Objects were grouped into two categories, "good" and "marginal." Marginal objects are 
questionable cases due to low S/N, bad sky subtraction, and most often, failure of A4959 to 
reproduce the intensity ratio of the components in A5007. Our objects are listed in [T] and 
their radio properties (see discussion below) are given in [2j 

Generally, the double-peaked profile was confirmed in lines such as [Ne ill] A3869 and 
the narrow H/3 component when the S/N was adequate. For the [S n] AA6716,6730, the 
typical [O ill] double peak splitting of ~ 400 or 500 km s" 1 is such that the A6730 of the 
lower redshift component falls on top of the A6716 line of the higher redshift component. 
This leads to a distinctive triple-peaked structure that was often seen when [S n] was within 
the spectral range. For [O n] AA3726, 3729, the profile rarely showed a double peak. The 
separation of the A3726 and A3729 lines is 225 km s -1 . For a typical [O ill] velocity split 
of ~ 400 km s -1 , the A3726 line of the higher redshift component fills in the gap between 
the A3726 and A3729 lines of the lower redshift component. This often gives a single peaked 
profile. To verify this, we computed a simulated [O il] blend for a typical candidate quasar 
using the relative intensity and velocity separation of the double-peaked [O ill] profile. The 
simulated [O il] profile does not show the double peak. (In a few cases, the splitting is 
seen in [O il]. The most dramatic case is in J0126+14, where the [O il] is clearly resolved 
into two individual lines with the same 580 km s -1 velocity splitting as in [O ill]). Because 
[O il] is almost never resolved into two components, [O ill] is split in all of our objects, and 
[S il] often exhibits a triple-peaked structure, a distinctive 1-2-3 pattern is often seen in our 
objects if [S il] is within the spectral range: one peak on [O il], two peaks on [O ill], and 
three peaks on [S il], as seen in Figured 

To confirm that both velocity components are typically present in [O il], we measured 
the centroid wavelength of the [O il] feature and compared this redshift with the redshift of 
the blue and red components of the double-peaked [O ill] profile. For a control sample of sev- 
eral normal QSOs with single-peaked [O ill], the [O il] and [O ill] redshifts agreed to within 



15 km s 1 , assuming a typical doublet ratio 7(A3726)//(A3729) = 0.83 (ISalviander et al. 



20071 ). Very few of our double-peaked [O ill] objects were consistent with an [O il] doublet 
being present solely at either the red or blue [O ill] redshift. For the double peaked [O ill] 
objects, the [O il] centroid was uniformly distributed between the red and blue [O ill] red- 
shifts, consistent with the fact that the primary (stronger) [O ill] peak can be either the red 
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or blue peak with similar frequency. The [O il] centroid redshift was on average 37% of the 
way from the primary [O ill] redshift to the secondary [O ill] redshift, implying that the 
stronger [O Hi] c omponent corresponds to the stronger [O II] component in most objects. 



Zhou et al.l (l2004f ) note the absence of a double peak for [O il] in most double-peaked [O ill] 
AGN, and suggest that this can be reconciled with a binary nature for these objects if the 
[O il] region is spatially more extended than for [O ill]. Our results suggest that [O il] does 
typically show both [O ill] velocity components, and no such explanation is necessary. 

Many of our objects, including a number of the AGN 1, show stellar absorption lines 
from the host galaxy in their spectrum. This affords an opportunity to exa mine the velocity 



of t he [O ill] compone nts relative to the velocity of the host, as done by iLiu et al.l (120091 ) 



and IWang et al.l (120091 ) for their AGN 2 with double-peaked [O ill] . In particular, Wang et 
al. offer a possible "virial" test of binarity involving the relative fluxes and velocities of the 
[O ill] components. In order to estimate fluxes of the [O ill] components, and to have a 
measure of the velocities independent of the cursor measurements described above, we carried 
out double Gaussian fits to the spectra of our double-peaked objects. The procedure assumed 
two Gaussians for the two components of the A5007 and A4959 lines, with adjustable velocity, 
FWHM, and flux, and the 3:1 intensity ratio required by atomic physics. The continuum 
in the vicinity of [O III] was fit with a second-order polynomial. The fit was optimized 
in a least-squares fashion. This procedure achieved a good fit to the total observed line 
profile in 106 of the 148 objects. For some of these objects, however, a visual inspection 
suggests that there could be substantially different choices for the amount of flux in the 
two components. For 85 of the 106 successful fits, we judged that the allocation of flux 
to the two components was unambiguous; these values are given in in Table [3j We found 
reasonable agreement between the velocities measured with SPLOT and the Gaussian fits. 
For example, the difference between the SPLOT and Gaussian redshifts of the red [O ill] 
component had a mean of 16 km s _1 and an r.m.s. scatter of 50 km s -1 in the rest frame. In 
most of the cases with a larger discrepancy, inspection of the fits suggested that the SPLOT 
measurements were more reliable, perhaps because the [O ill] line is typically more sharply 
peaked than a Gaussian. We have used the SPLOT velocities below; use of the Gaussian 
velocities would not alter our conclusions. 

Host galaxy velocities were determined by cross-correlating a template galaxy spectrum 
with the observed spectrum with the aid of the IRAF routine FXCOR used in interactive 
mode. All 148 of the double [O ill] spectra were visually examined in the region containing 
the MgB and NaD lines. If these features were visually detected (73 objects), the spectrum 
was tagged and converted to a linear wavelength scale spanning 3600-8500 A with dispersion 
of 2.0 A/pixel. The cross-correlation template was generated with the same wavelength scale 
and dispersion, using an early type galaxy SDSS spectral template. FXCOR was used to 
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cross-correlate each galaxy spectrum with the template. A cubic spline was used to model 
and subtract out the mean continuum shape and a cosine bell was used to apodize 1% of 
each ID spectrum at both ends. Interactive options were used to cross-correlate the spectra 
only around the MgB-NaD region of the spectrum. If a well defined cross-correlation peak 
was not visible, the cross-correlation output was rejected. This procedure was successful for 
67 objects, of which 46 have [O ill] component fluxes in Table [31 Table [1] gives the host 
galaxy redshift. The host redshifts have an average uncertainty of 45 km s _1 in the rest 
frame, with a maximum uncertainty of 95 km s" 1 . 



3. Double-Peaked [O ill] Objects 

We have a total of 148 double-peaked [O ill] objects, listed in Table [TJ Objects with 
an asterisk exhibit either binarity or signs of interaction in the Sloan image. Also given are 
the SDSS redshift, the velocity separation of the [O ill] peaks, the continuum luminosity 
vL u (5100) , and the AGN type (1 if broad emission lines are visible). The "Quality" column 
denotes whether the object is good, g, or marginal, m. There are 78 good and 70 marginal 
objects. While marginal objects do have a listed velocity splitting, these measurements are 
uncertain in some cases. Figure [1] gives spectra for three representative objects. 

Figure El gives the SDSS multi-color images of the eight visually interesting objects. 
Two of these, J1316+17 and J1441+09, are included because they exhibit numerous smaller 
companions within a close radius. Three of the images appear to consist of two or more 
individual galaxies at very close proximity: J 1245+37 and J1307+46 appear to show pairs 
of galaxies centered on the SDSS fiber, and J1157+08, shows three galaxies in a linear 
arrangement. 

J0942+12 exhibits large tidal tails in the image, suggesting recent interaction. However, 
the nucleus itself does not appear binary. 

J1001+28 may be a major merger between a large spiral and a spiral or elliptical. While 
the fiber width is not enough to encompass the NLRs of both objects, the galaxy is clearly 
undergoing interaction. 

J1516+18 (not shown) has 1.3 Jy of radio flux, an order of magnitude larger than any 
other object in our sample. The object is a good double [O ill], although the velocity split 
which is clear in [O ill] is not seen in other lines. 



J1517+33 is discussed at length in lRosario et al.l (120 101 ). It has emission line knots that 
straddle the stellar nucleus in velocity and position. The double peaked profile appears to 
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be the result of a bipolar radio jet. 



Statistical Results 



We did not impose an a priori S/N cutoff, in order to avoid missing noisy but interesting 
objects with strong [O ill] lines. However, some of the SDSS spectra did not have sufficient 
S/N to show the typical [O ill] double peaks displayed by our objects. Therefore, for statis- 
tical analysis, we must estimate what fraction of the full quasar sample had sufficient quality 
to be considered the parent sample of our double-peaked objects. We find by inspection of 
a substantial subset of the 21,592 objects that only 40% ± 10% of spectra are viable in the 
sense that the double peaks of A4959, as present in a typical "good" object, would be clearly 
distinguishable from the noise. While subjective, this estimate is reasonably sound because 
of the strong contrast between the S/N of the better and worse spectra. Moreover, there are 
larger uncertainties in the analysis below, such as corrections for binaries with line-of-sight 
velocities too small to give a double-peaked [O ill] profile (see below). 

An issue for statistical analysis is the presence of some AGN 2 in the SDSS quasar 
download, numbering about 2.1% of the full sample. Interestingly, 42% of our objects are 
narrow-line objects, so that the Type 2 objects are twenty times overrepresented among our 
double [O ill] objects. We believe this to be a selection effect. SDSS defines a quasar as 
having one or more emission lines with FWHM of at least 1000 km s" 1 , intended to eliminate 
narrow-line objects. Evidently, the pipeline failure responsible for the presence of some AGN 
2 in the parent sample is influenced by a double-peaked profile. Because of this dramatic 
selection effect, we omit the Type 2 objects from our statistical analysis. 

After subtracting 2.1% for the Type 2 objects, we have a parent sample of 8,452 viable 
Type 1 spectra. Our 86 Type 1 double [O ill] objects amount to 1.0% of the parent sample; 
42 of these (0.50% of the parent sample) are "good". Our incidence of 1.0% double-p eaked 
AGN 1 is similar to the incidence of 1.1% for AGN 2 in the study of iLiu et all (120091 ). 



Our double-peaked AGN were cross-referenced with the FIRST radio catalogue ((Becker et al 



19951 1 . Only four of our objects were outside the footprint of the FIRST survey. FIRST 
detects 20 cm radio flux for 27% of our Type 1 objects (detection limit 1 mJy), as com- 
pared to only 9% FIRST d etections in the overall SDSS quasar catalogue in our redshift 
range ((Schneider et al.l 120071 ). In other words, radio detection is three times overrepresented 
among our double-peaked [O ill] objects. This suggests that radio jet interaction is often 
the c ause of the double-peaked [O ill] profiles. Su ch interactions in the NLR are well known 



(e.g. : 



Whittle fc Wilson 



20041 ; IWhittle et al.ll2005l and references therein) . One of our double 
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[O ill] objects, J1517+33, appears to involve such an interaction. It is a narrow-line object 
with a FIRST radio flux of 107 mJy and two bright optical regions near its nucleus in the 
Sloan image. Subsequent observations of this object with the VLA have shown a double- 
lobed radio structure at the same or i entati on and position as the bright optical regions. This 



object is discussed in iRosario et al.l (120101 ) . 



We give the radio fluxes, luminosities, and the iKellermann et al.l (I1989|) radio loudness 



param eter in Table [2J We adopt the radio-loudness conventions used in 



Kellermann et al. 



( 19891 ). The ratio R = F r /F measures the ratio of 6 cm radio flux density to 4400 A optical 
flux density. If 0.1 < R < 1, we call an object radio quiet, 1 < R < 10 is radio intermediate, 
and 10 < R < 100 is radio loud. If an object is not detected in the FIRST sample, we refer 
to it as radio undetected. It should be noted that the lack of a FIRST detection does not 
indicate a total absence of radio flux, or of radio jets. High redshift objects with no FIRST 
detection may have marginal flux which could have been detected at lower redshifts, and 
orientation of the radio source can affect detected flux. Therefore some objects with radio 
jets may fall below the FIRST detection limit. 

Additionally, Type 2 spectra are dominated by the light of the host galaxy, which 
clouds the meaning of the radio-to-optical flux. Consequently, the R ratios given in Table [2] 
for Type 2 objects should be taken with caution. Because our study focuses on broad-line 
objects and we omit the AGN 2 from our statistical interpretation, this issue does not affect 
our conclusions. 

Because of the likelihood of radio jets causing the double-peaked signature, we consider 
our best candidates for true binaries to be those with no FIRST radio detections. Discarding 
9% of the parent sample to allow for FIRST-detected objects, our parent sample consists of 
7,692 quasars. Our double-peaked sample consists of 66 radio undetected broad-line objects, 
those we consider to be the best prospects. This is 0.9% of the total number of viable SDSS 
radio undetected broad-line quasars. If we use only the "good" objects, 29 remain, which is 
0.38% of the parent sample. 

Figure H] shows the relationship between the redshift and the velocity split for our 
double-peaked AGN. Objects at greater distance ar e typically at higher l uminosity, and 



more luminous objects tend to have wider [O ill] lines (ISalviander et al.ll2007l and references 
therein). Accordingly, a larger velocity difference would be necessary for a resolved double- 
peaked profile. It would make sense, therefore, that as redshift increases, observed velocity 
splits should also increase. However, we find that the opposite is true. It appears that at high 
redshifts, the typical velocity splitting actually decreases, although the statistical sample is 
small. 
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When the redshift plot is made separately for AGN 1 and AGN 2, a distinctive gap is 
seen among the AGN 2 between redshifts 0.3 and 0.5. In a control sample of the SDSS AGN 
2 in our quasar download, we also see a gap in that redshift range. We assume therefore 
that the gap results from the failure in the SDSS pipeline that allowed some AGN 2 into the 
quasar sample in the first place. 

For our entire set of double-peaked objects, the redshifts range from z = 0.100 to 0.686 
with a mean of 0.33. For our "good" AGN 1, the average EW is 86 A. For the Type 1 objects 
with host galaxy redshifts, the midpoint in velocity between the red and blue components 
averages +33 ±29 km s _1 (standard error of mean), with an r.m.s. scatter of 108 km s -1 . For 
our AGN 1 with [O ill] component fits, the mean value of the Gaussian width parameter for 
the red component is cr r = 143 ± 15 km s _1 , with an r.m.s. scatter of 57 km s _1 ; for the blue 
component, the mean is = 161 ± 18 km s _1 , with a scatter of 66 km s _1 . For a Gaussian, 
the FWHM is 2.35a. For this group of objects, the mean separation is A V = 364 km s -1 



which is 1.02 times the mean FWHM of either component. For comparison, iLiu et al.l (120091 ) 



found a typical AV larger by a factor ~ 1.5 than the FWHM of each component. 



5. Discussion 

5.1. Reliability of Double [O ill] 

The "good" objects are definite cases of double-peaked [O ill] profiles. The split is 
confirmed with good S/N and velocity agreement in A4959, and in other narrow lines when 
there is adequate S/N. For the marginal cases, the confirmation is less compelling in A4959, 
and the S/N may be inadequate in other lines. In order to estimate the frequency with 
which noise might give the appearance of a double-peaked [O ill] profile , we examined a 



set of simulated single-peaked spectra for 1,773 SDSS quasars computed by lSalviander et al. 



(120071 ) as a test of their procedures. These objects are typical of our double-peaked sample in 
terms of [O ill] strength and continuum S/N. The simulated spectra reproduce the observed 
emission lines fo r each quasar, using a G auss- Hermite fit to the lines in the original spectrum 



as described by ISalviander et al.l (120071 ). along with a power-law continuum approximating 
the observed continuum in the [O ill] region. Gaussian random noise was added pixel- 
by-pixel using the noise amplitude given in HDU3 of the SDSS spectrum. For the present 
purposes, the simulated spectra were visually inspected in the same manner as our actual 
SDSS parent sample. In no case did a double-peaked A5007 profile appear as a result of 
noise that would have met our standards even for a "marginal" object. We conclude that all 
of the "good" and most of the "marginal" objects are in fact genuine double-peaked [O ill] 
AGN. 



- 10 - 



5.2. Binary AGN? 

Our search was motivated by the prospect of discovering binary AGN. How many of our 
objects are actually binaries? The Type 1 objects have broad emission lines and are clearly 
AGN. The Type 2 objects have large equivalent widths of [O ill], a typical [O lll]/H/3 
intensity ratio ~ 10, characteristic of AGN, and sometimes show He II A4686 or [Ne v]. This 
indicates that most are indeed AGN. As noted above, alternative interpretations of double- 
peaked [O ill] profiles include jet interactions and rotation. One of our Type 2 objects, 
J15 17 +33, is spatially res olved but is clearly an example of a radio jet interaction with the 



NLR (IRosario et al.ll2010l ). Another of our objects, J1129+60, was observed with the VLA 
by Rosario et al. It also has kpc scale radio structure and may involve jet interaction. The 
majority of the radio undetected objects show double-peaked line profiles of a type that could 
reasonably occur for a binary AGN. One AGN 1, J1307+46, shows two spatially resolved 
optical sources separated by 3 arcsec (15 kpc). The others are unresolved in the SDSS images 
and may be candidates for binary AGN at spacings less than about 5 kpc, depending on 
redshift. High resolution imaging in the optical, radio, and X-ray is needed to help determine 
which are actually binary AGN. Pending such studies, can other arguments offer guidance 
as to which objects are binaries or disturbed NLRs, or at least indicate the proportions of 
such objects in our sample? 



5.2.1. Velocity separation 



The velocity separation of the double peaks is plausible for binary AGN in most cases. 
For a circular orbit in the gravitatio nal potential of an isotherm al sphere of velocity dispersion 
cr*, the orbital velocity is 2 1 ' 2 a* (IBinney fc Tremaind 119 88). For an eccentric orbit, the 
pericentric velocity could be several times ex*. Our "good" radio- undetected Type 1 objects 
have a mean value log Mbh = 8.13 solar masses, based on the broad line width and continuum 



luminosity using equation 2 of lShields et all (120031 ). (This refers to a subset of 13 objects with 
the best determined Mbh-) This corresponds to a host bu lge velocity dispersion l og a* = 2.30 
km s _1 based on the Mbh — c* relationship as given by iTremaine et al.l (120021 ). The mean 
[O ill] velocity split is log Aw = 2.52 in km s _1 ; Av/a* averages 1.7 with a range of 1.0 to 
6. A precise prediction of Ay/a* for binary AGN with double [O ill] is difficult because 
of unknown orbital parameters and a selection for binaries whose relative velocity vector is 
oriented parallel to the line of sight. In any case, the largest separations that we see, up to 
1400 km s _1 , seem rather high for mergers and may represent cases of a bipolar jet. 
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5.2.2. Line intensity ratios 

The utility of emission line ratios to distinguish binaries from disturbed NLRs is unclear. 
Jet interaction sources may involve gas dynamically affected by the jets but stil l photoionized 



by th e AGN continuum and showing normal photoionization line ratios (e.g.. IWhittle et al. 



20051 ). Binarity could conceivably affect the NLR line ratios, particularly from spacings 
< 1 kpc approaching the size of the NLR itself. For our radio-undetected "good" Type 1 
objects, average values are EW(A5007) = 72 A and F(A3727)/F(A5007) = 0.24, where A3727 
refers to the combined intensity of the [O n] doublet. These values are not significantly 
different from a control sample of comparable non-double [O I II] SPSS quasars, an d the 



[O n]/[0 ill] ratio is typical of power-law photoionized objects ([Baldwin et al.lll98ll ). The 
[O i] A6300 line is measureable in only a fraction of our objects. For those objects, the 
average value of 7(A6300)//(A5007) is 0.11, similar to a comparable control sample; and these 
ob jects lie in the power -law photoionized region of the [O n]/[0 ill]- [O i]/[0 ll]diagram 



of Baldwin et al.l ( 119811 ). For our radio-detected "good" Type 1 objects, average values are 
EW(A5007) = 128 A and F(A3727)/F(A5007) = 0.17. The greater strength of the narrow 
emission lines in radio loud AGN is well known. The mean values for [O n]/[0 ill] differ 
by less than 2 a between the radio-detected and undetected samples. Given the impact of 
radio emission on the narrow lines for normal (non-double peaked) AGN, comparisons of 
line intensities for radio-detected and non-detected objects may not give clear guidance as 
to which double-peaked objects are binaries or jet-interactions. 

Differential reddening of the red and blue components of the double-peaked narrow lines 
may offer a diagnostic (M. Whittle, private communication). Such a reddening differential 
is expected for bipolar jets but not for binary AGN or a rotating disk geometry. If there is 
distributed reddening in the NLR, the blueshifted component, being on the near side of the 
AGN, may show less reddening than the redshifted component. Allowing for noise, blending, 
and redshift, only a handful of our objects would allow a reliable assessment of the Ha/H/3 
intensity ratio for the two components separately, too few for a meaningful test. For our 
objects with measured [O ill] component fluxes, the red/blue ratio is F r /F h = 0.90 ± 0.10 
for the AGN 1 and 1.30 ± 0.20 for the AGN 2 (standard error of the mean). These values 
give no clear indication of a systematic difference in extinction, but the uncertainty is large, 
reflecting the scatter among individual objects in F r /F^. 



In a study of narrow line AGN in SDSS with double-peaked [O ill] profiles, iLiu et al. 



( 120091 1 find that the intensity ratio of [O ill] to H/3 is typically rather similar in the red and 
blue velocity components. This is suggestive of a common ionizing source and may favor jets 
or disks over binary AGN for most objects, at least in their sample. 



- 12 - 



5.2.3. Virial test 



In an independent study of AGN 2 with double-peaked profiles, I Wang et al.l (120091 ) 
examine the relative velocity offset of the [O ill] red and blue components (v T /vb) in relation 
to their relative intensities (F r /F^). They find that the brighter component tends to be closer 
to the velocity of the host galaxy. They argue that this is consistent with orbital motion of a 
binary galactic nucleus, assuming that the brighter emission-line component typically comes 
from the more massive black hole (and associated stellar cluster) that would normally have 
the lower orbital velocity, F v /F\> oc (t> r /t>b)~ 1 - 

We have examined our own data set for the trend found by Wang et al. We considered 
separately the AGN 1 and AGN 2, using flux determinations for the two [O ill] line compo- 
nents as given in Table 3. The results are shown in Figure [5j There is little indication of the 
claimed virial trend for our AGN 1 sample. Unfortunately, the number of useful AGN 1 is 
only 13. This includes 7 objects with v v /v\ } > 1 and 6 objects with v r /v^ < 1. For these AGN 
1 and AGN 2, respectively, the average values are f r A'b = (2.33, 0.65) and F r /F h = 0.92, 1.04. 
These values are consistent with no dependence of velocity offset on flux ratio at the 1 a 
level, but they differ by 4 a from the reciprocal relation predicted by the virial argument of 
Wang et al. For our AGN 2, Figure [5] suggests some inverse relationship between the velocity 
and flux ratios, considerably weaker than (v T /v^)~ 1 . However, the sample is small (23 useful 
objects), and our AGN 2 sample is subject to unknown biases, as discussed above. 



5.2.4. Luminosity dependence 



One of our ob j ects, SPSS J131642. 90+175332. 5, was the subject of a detailed study 
by Ku &: Komossal ( 120091 ) . These authors consider alternative explanations of the double- 
peaked lines, including a binary AGN, biconical outflow, jet-cloud interaction, and other 
complexities in the NLR geometry. They suggest that binaries might be more common 
in quasars than in Seyfert galaxies, on the assumption that major mergers are more often 
involved in fueling quasars. We have examined our sample for such a trend. For luminosity 
bins log uL u (U00 A) in the ranges (<43.5, 43.5 - 44.0, 44.0 - 44.5, 44.5 - 45.0, > 45.0), the 
number of AGN l's in our sample is (4, 20, 34, 24, 4). These counts represent (0.35, 0.30, 
0.39, 0.74, 1.13) percent of the number of objects in the parent sample in the corresponding 
l uminosity bins . For reference, the Seyfert /quasar boundary is Mb = —21.5 + 5 log h 
( Peterson! Il997 ). corresponding to log z/L^(4400 A) = 44.44. The incidence of double peaks 
increases roughly a factor of two across our range of luminosities. This resembles the factor 
of two in crease in the inc idence of double-peaked [O ill] with increasing [O ill] luminosity 
found by lLiu et al.l (120091 ) in their study of AGN 2. This could be indicative of a substantial 
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number of binaries among our sample. However, the origin of bipolar jets in AGN is not 
fully understood, and jets affecting the NLR could be more common in higher luminosity 
AGN. 



5. 2. 5. AGN 1 versus AGN 2 

We suggest here another statistical argument involving the velocity splitting that may 
give an indication of the nature of the double-peaked narrow line AGN. This involves a 
comparison of the typical velocity separation for Ty pe 1 (broad lines) and Type 2 AGN 



(narrow lines only). In the unified model of AGN (lUrry fc Padovanil 119951 ) . AGN 1 are 
viewed relatively close to the disk axis, whereas AGN 2 are viewed closer to the disk plane. 
In AGN 2, the "dusty torus" obscures the central continuum source and BLR. Thus, the 
typical velocity projection onto the line of sight may differ between AGN 1 and AGN 2. 
(1) In the bipolar jet picture, the observed splitting is a function of the observer's location 
relative to the jet axis, V = Vocos# where Vq is the true velocity separation. In the unified 
model, on average, AGN 1 should have larger observed velocity separations than AGN 2. 
Consider a simple model in which all objects have the same Vo, and in which all objects 
have a torus opening angle of 45 degrees. Let it?i, w 2 represent the average observed velocity 
splitting for AGN 1 and 2 respectively. Then for AGN 1, W\ — Vq(cos8), where the average 
is over polar angles of to 45 degrees; and for AGN 2, the range of 9 is from 45 to 90 degrees. 
This gives (cos0) = 0.85 for AGN 1 and 0.35 for AGN 2, so that wi/w 2 = 2.4. However, this 
exaggerates the effect, because there is a minimum splitting for the double-peaked profile 
to be resolved. Our results suggest a value for this of about 200 km s -1 . For a rough 
estimate of this correction, let us take Vo = 500 km s _1 , based on the mean splitting for 
our AGN 1. Thus, AGN with cos# < 0.4 will not show double peaks. The range of polar 
angles for AGN 2 to show double peaks is then 45 to 66 degrees, giving (cos 6) = 0.55 and 
W1/1V2 = 1.54. (2) In a rotating disk model with circular velocity V rot , the observed velocity 
splitting is V = Vosin^, where V = 2V TOt . By a similar procedure to that above, we find 
u>i/u>2 = 0.53 for V min = or W\/w 2 = 0.63 using V min = 200 km s _1 . Thus, the sense of 
the difference between AGN 1 and 2 is reversed from the jet case, but the magnitude of the 
effect is similar. (3) For binaries, the prediction is less clear. Observationally, only one of 
the two AGN need be a Type 1 for broad lines to be observed. If the dusty tori of both 
objects tend to be coplanar with the orbit, then the situation resembles case (2) above for 
a rotating disk, and AGN 1 might show smaller velocity separations than AGN 2. If there 
is little correlation between the orbital plane and the dusty tori, then w\ and w 2 should be 
similar. 
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Our observed double-peaked [O ill] objects have w\ = 434 ± 25 km s" 1 for AGN 1 
and 480 ± 22 km s -1 for AGN 2 (standard error of the mean). The AGN 1 and 2 samples 
have closely similar mean luminosities in [O ill], consistent with similar underlying AGN 
luminosities. The observed ratio of W\jwi = 0.90 ± 0.07 disagrees with a predominance of 
either bipolar jets or disks, if the above estimates are roughly correct. This result, by itself, 
might suggest either similar numbers of disks and jets, or a predominance of binaries. We 
refrain from drawing conclusions, however, because of possible biases related to the presence 
of Type 2 AGN in our sample. 



The double-peaked objects of iLiu et all ( 120091 ). all AGN 2, have W2 = 368 ± 8 km s _1 , 
close to our value. Our sample has a mean redshift of z — 0.37 and 0.27 for AGN 1 and 2, 
respectively. For Liu et al. the mean redshift is 0.16, and only 0.088 for Wang et al. For a 
subset of our lower redshift type 1 AGN having an average redshift of 0.158, similar to that 
of Liu et al., we find that our average velocity splitting is w\ = 425 ± 58 km s" 1 , within 1 a 
of Liu et al.'s sample, and that w\jw-i = 1.1, again consistent with binaries. If we compare 
Li u et al.'s AGN 2 to our full sample of Type 1 objects, we find w±/w2 = 1.2. The objects 
of IWang et all (120091 ) have u>2 = 344 km s -1 , but their sample has a much smaller average 
redshift of 0.088, and their mean [O ill] luminosity 10 41 - 30 erg s _1 is an order-of-magnitude 
smaller than for our sample. 



5.3. Line-of-sight correction 



The width of the narrow emission lines correlate s with a* empirically, with log ofo ttt]/o", 
0.0 ± 0.15, where a [Q m] = FWHM(A5007)/2.35 (kelson fc Whittle! Il996t bonning et al. 



20051 ). The profile will appear double-peaked if the velocity separation is greater than 
~ o"[o in]- Given that the orbital velocity of a binary should be of this order (see above), one 
might expect that some binaries will show a double-peaked [O ill] profile and others will not. 
Moreover, a substantial fraction of binaries may be at a low velocity phase of an eccentric 
orbit, and many will have their relative velocity vector close to the plane of the sky. The 
true number of binaries (or otherwise distur bed NLRs) may b e several times larger than the 
number of double-peaked objects observed. IZhou et all (120041 ) estimate that this correction 
may be an order-of-magnitude. Thus, if a large fraction of our radio-undetected Type 1 ob- 
jects are binaries, then the implied true incidence of binaries is as large as ~ 10%. However, 
such a large fraction of binaries appears to conflict with the statistics of spatially resolved 
AGN (see below). This in turn suggests that only a small fraction of the double-opeaked 
[O ill] objects are in fact binary AGN. 
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5.4. Fueling of Binary AGN 



The role of mergers in fueling AGN activity is an important topic of study. While 
seemingly undisturbed AGN are observed, high resolution imaging of quasa rs shows close 



companions and and feat ures of gravitational interaction in many cases (e.g., iBahcall et al. 



19971 ; iBennert et al.ll2008l ). Simulations of galactic collisions with gas show that tidal torques 
lead to concentrations of ga s in each nucleus afte r the first close encounter, and in the nucleus 
of the final merger product (IHopkins et al.ll2006l . and references therein). It seems likely that 
a large fraction of luminous AGN are triggered by mergers. This suggests that binary AGN 
would be a common occurrence. 

As noted above, the incidence of known close binaries at 1 or 2 kpc separation is only 
~ 10~ 3 , and it has been ar gued that this is lower th an expected if both black holes in a merger 
are fueled inde p enden tly fjJunkkarinen et al.ll200ll ). A similar conclusion is also reached by 
Foreman et al.l (120091 ). in a study of optically resolved binaries in the 10 to 100 kpc range. 
In our case, a typical redshift is z ~ 0.3, and a 1 arcsec limit on our unresolved objects 
corresponds to ~ 5 kpc. What might one expect for the number of binaries in this range 
of separation? For an orbital velocity of 300 km s" 1 and radius of 3 kpc, the orbital period 
is t or b ~ 10 7 ' 8 yr. Dynamical friction causes the orbit to decay o n a timescle of a few 
times t orh ; for typical parameters for our objects, t d{ m 10 8,9 yr (IJunkkarinen et al.l 120011 ). 
The lifetime of an AGN outburst is often estimated to be the Salpeter growth timescale 
£agn ~ 10 7 ' 6 yr. Thus, the duration of an AGN outburst could be a substantial fraction of 
tdf- The fact that the galaxy is observed as an AGN implies that at least one black hole is 
currently fueled. If fueling of each hole occurs independently, then the probability that a 
given AGN will have both holes fueled as a binary AGN is p^n ~ ^AGNAdf ~ 10~ L3 . This 
assumes that most AGN are triggered by mergers. 



This theoretical estimate is considerably larger than our observed ~ 1% incidence of 
double [O ill] objects but similar to the above 10% estimate after correction for unresolved 
double [O ill] objects. It is doubtful that the true incidence of binary AGN is so high, 
however. The number of resolved binary QSOs at the 2 kpc scale and the 30 kpc scale is 
~ 10~ 3 in both cases , based on objects at typical redshifts of 1 to 3 (IJunkkarinen et al.ll2001 



Foreman et al.ll2009i ). If these indications of a low binary rate carry over to our redshift and 
separation range, this may imply that most of the double [O ill] objects are in fact single 
AGN with a disturbed NLR. A low number of optically visible binary AGN could result 
from a delay in fueling of the black holes until the orbital separation has decayed below the 
kpc scale, or obscuration of m ost binary AGN at kpc separations by massive amounts of 
dusty gas (IHopkins et al.ll2006l ). In this context, it is interesting that ul traluminous infrared 



galaxies (ULIRGs) often show nuclear spacings of one or several kpc (ISanders &: Mirabel 
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1996). 



6. Conclusion 



We have found that approximately 1% of broad-line SDSS quasars at 0.1 < z < 0.7 have 
double-peaked [O ill] line profiles, and other narrow lines show similar velocity splittings 
when the data quality is sufficient. In most cases, the magnitude of the velocities is consistent 
with orbital motion in binary AGN with distinct NLRs. However, many of the objects, 
including those with FIRST radio detections, are likely disturbed NLRs involving bipolar 
jets or other complexities. We have discussed a number of potential statistical tests of 
which mechanism accounts for most of the double-peaked objects, but available data are not 
definitive. 

AGN with double-peaked narrow lines are interesting individual targets for detailed 
study. In addition, they offer a potential way to assess the probability of simultaneous 
fueling of both black holes in a galactic merger. In order to realize this potential, further 
study of these objects, including high resolution imaging in the radio, optical, and X-ray, is 
needed to determine which of them are actually binary AGN. 
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Fig. 1. — Spectra in the [O ill] and H/3 region for three examples of objects with double- 
peaked [O ill] profiles from Table 1. Figure 3 gives the SDSS image of J1307+46. In our 
full sample, the red and blue peaks are stronger with roughly equal frequency. See text for 
discussion. 
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Fig. 2. — The distinctive 1-2-3 pattern typical of many of our objects. The leftmost panels 
show the [O n] lines, nearly always single-peaked. The center panels display the double- 
peaked [O ill] lines. The last panels show the triple-peaked structure of the [S n] line, 
discussed in the text. 
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J1316+17 J1245+37 J 1441 +09 J 151 7+33 




J 1307+46 J 1001 +28 J0942+12 J1 157+08 



Fig. 3. — Eight SDSS images are displayed, corresponding to the eight objects in Table 1 with 
asterisks following their names. The scale bar in the upper left of the pictures corresponds 
to 5 arcseconds. The images exhibit objects with two or more obvious components, double 
nuclei, or signs of recent interaction. See text for discussion. 



-22 - 



1500 



oo 1000 
a 



Oh 

£ 500 
o 

u 

> 



1500 



1000 



• 


Radio Undetected 


▲ 


Radio Intermediate 


■ 


Radio Loud 


• 



A A - • • 



• ■. • 



_^ I L 



500 







0.1 0.2 0.3 0.4 



0.5 



0.6 



I 9 



• • • 



•.A . 



j_ 



j_ 



0.1 0.2 0.3 0.4 0.5 0.6 

Redshift 



Fig. 4. — Velocity splitting of the type 1 (upper panel) and type 2 double- [O ill] objects as a 
function of redshift. There appears to be a slight trend toward fewer large splittings at larger 
redshifts. The graph is color-coded for radio-loudness, where this parameter is determined 
by the Kellerman R ratio. The gap at z ~ 0.3 to 0.5 for the AGN 2 is an artifact. Note that 
radio detection by FIRST is redshift dependent. See text for discussion. 
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Fig. 5. — Flux ratio of the red and blue [O ill] components as a function of the ratio of 
their velocity offsets relative to the host galaxy redshift (see text). Upper and lower panels 
show AGN 1 and AGN 2, respectively. 
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Table 1. Double-Peaked AGN 



SDSS Name 


2SDSS 


LI 

(km/s) 


log vL v 

(erg/s) 


Spectral 
Type 


Quality 






•^host 


J011802.94-082647.2 


0.137 


350 


43.69 


2 


m 





.13771 


0. 


13713 


J012613.31+142013.4 


0.573 


580 


44.35 


2 


g 





.57338 






J020011. 52-093126.1 


0.361 


200 


44.49 


1 


m 





.36077 






J072554.42+374436.9 


0.634 


250 


44.81 


1 


m 





.63500 






J074129.66+392835.9 


0.210 


480 


43.63 


2 


m 





.21059 


0. 


20982 


J080315.67+483603.1 


0.635 


280 


44.86 


1 


m 





.63576 






J080841. 21+481351.9 


0.123 


340 


43.38 


2 


m 





.12388 


0. 


12359 


J081507.41+430427.0 


0.510 


500 


44.09 


1 


g 





.51014 






J081542. 53+063522.9 


0.244 


280 


44.12 


1 


m 





.24428 






J082357.80+391630.9 


0.166 


700 


43.67 


1 


m 





.16774 


0. 


16630 


J082857.99+074255.7 


0.554 


300 


44.07 


2 


g 





.55475 






J084049.46+272704.7 


0.136 


580 


43.32 


2 


m 





.13706 


0. 


13559 


J090246.93+012028.2 


0.513 


240 


43.88 


2 


m 





.51408 






J090615.92+121845.6 


0.644 


200 


45.23 


1 


g 





.64410 






J090947.85+312443.6 


0.264 


1050 


44.41 


1 


g 





.26449 






J090958. 32+085542.2 


0.158 


240 


43.20 


2 


m 





.15830 


0. 


15869 


J091110. 19+140632.6 


0.387 


350 


44.87 


1 


g 





.38744 






J091459. 25+210219.6 


0.133 


540 


43.32 


2 


g 





.13418 


0. 


13335 


J091649.41+000031.5 


0.222 


310 


43.47 


1 


m 





.22307 


0. 


22223 


J091654.09+521723.0 


0.219 


315 


43.70 


1 


m 





.21998 






J092152.46+515348.0 


0.587 


440 


44.06 


2 


m 





.58806 






J092455. 24+051052.0 


0.149 


410 


43.85 


2 


m 





.15005 


0. 


15042 


J094100. 81+143614.4 


0.383 


870 


44.20 


2 


m 





.38493 






J094144.83+575123.6 


0.159 


300 


43.69 


1 


g 





.15924 


0. 


15866 


J094205. 83+125433.6* 


0.154 


680 


43.95 


2 


g 





.15443 


0. 


15420 


J094236.68+192541.1 


0.540 


480 


45.23 


1 


g 





.54249 






J095207.62+255257.2 


0.339 


475 


44.35 


1 


g 





.33986 


0. 


33892 


J100145.30+283330.3* 


0.114 


550 


43.15 


2 


g 





.11543 


0. 


11550 


J100708.01+242039.0 


0.544 


250 


44.62 


1 


g 





.54497 






J101034.28+372514.7 


0.282 


1400 


43.98 


1 


g 





.28228 






J101241. 20+215556.0 


0.111 


210 


43.51 


1 


m 





.11126 


0. 


11069 


J102004.36+324342.5 


0.484 


340 


44.26 


1 


g 





.48484 






J102038. 74+212806.5 


0.137 


660 


42.53 


2 


m 





.13925 






J102045. 58+030306.9 


0.535 


260 


44.54 


1 


m 





.53540 






J102727.90+305902.4 


0.124 


360 


43.81 


2 


g 





.12488 


0. 


12451 


J103138.67+380651.7 


0.492 


730 


44.82 


1 


g 





.49211 
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Table 1 — Continued 



SDSS Name 


2SDSS 


V 

(km/s) 


log vL v 
(erg/s) 


Spectral 
Type 


Quality 






•^host 


J103326.92+284751.0 


0.591 


390 


44.66 


2 


g 





.59217 






J103752. 23+312500.3 


0.160 


900 


43.66 


2 


m 





16275 





.16020 


J104257.96+385347.2 


0.107 


440 


43.38 


2 


g 





40716 





.10661 


J105104.54+625159.3 


0.436 


360 


44.99 


1 


m 





.43734 






J111013. 20+053338.8 


0.152 


830 


43.46 


2 


m 





.15271 






J111710.07+333950.3 


0.128 


560 


43.69 


2 


m 





.12809 





.12786 


J112019.61+130320.0 


0.314 


600 


44.13 


1 


m 





.31486 






J112319. 21+302825.4 


0.522 


480 


43.50 


2 


g 





.52294 






J112507.33+023719.1 


0.260 


400 


43.92 


2 


m 





.26095 





.25989 


J112634.84+455935.7 


0.278 


300 


43.86 


2 


m 





.27884 





.27849 


J112939. 78+605742.6 


0.112 


550 


43.61 


2 


g 





.11276 





.11171 


J113020.99+022211.5 


0.241 


330 


44.19 


1 


m 
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Note. — AGN with double-peaked narrow emission lines from the Sloan Digital Sky Survey, in- 
cluding redshift as given by SDSS, rest-frame velocity separation of the two peaks of [O in] A5007, 
optical continuum luminosity from the SDSS spectrum, AGN type, quality of the double peak as 
described in Section 2, redshift of the red peak of A5007, and host galaxy redshift (when measur- 
able). Objects with an asterisk following their name are accompanied by an image in Figure [3l See 
text for discussion. 
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Table 2 — Continued 
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Note. - Radio properties of the double-peaked AGN that have detections in 
the FIRST survey. Columns give SDSS redshift, observed 20 cm flux density from 
FIRST, source specific lumi nosity at 6 cm rest wav elength (assuming L 
radio loudness parameter R ( iKellermann et al.lll989l ). and AGN type. 
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Table 3. A5007 Component Fluxes 
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Table 3 — Continued 
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Note. Fluxes of the red 

and blue components of the A5007 
emission line. Units of flux are 
10 -17 erg cm -2 s^ 1 . 



